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ABSTRACT 

We report spectral time series of the late O-type runaway supergiant HD 188209. Radial 
velocity variations of photospheric absorption lines with a possible quasi-period 6.4 
days have been detected in high-resolution echelle spectra. Night-to-night variations in 
the position and strength of the central emission reversal of the Ha profile occuring 
over ill-defined time-scales have been observed. The fundamental parameters of the 
star have been derived using state-of-the-art plane-parallel and unified non-LTE model 
atmospheres, these last including the mass-loss rate. The derived helium abundance is 
moderately enhanced with respect to solar, and the stellar masses are lower than those 
predicted by the evolutionary models. The binary nature of this star is not suggested 
either from Hipparcos photometry or from radial velocity curves. 
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1 INTRODUCTION 

Runaway O stars have been defined as a group by Blaauw 
(1961), who introduced the term runaway to describe the 
space motions of AE Aur and /j. Col. Blaauw (1961) has also 
suggested that such stars were ejected in the breakup of bi- 
nary systems in supernova explosions by their companions. 
In later evolutionary stages, the initial secondary appears 
as a most massive star and transfers matter to the compact 
companion (the initial primary) making the system appear 
as a massive X-ray binary (van den Heuvel 1976) . Given the 
possibility of the binary nature of runaway stars, it appears 
to be an important task to measure the radial velocity (RV) 
variations of the photospheric lines. Systematic searches for 
RV variations have been made in order to assess the bi- 
nary frequency of O stars (e.g. Garmany, Conti & Massey 
1980; Stone 1982, Gies 1987). In many CEises the amplitude 
of RV variations is quite large, and the additional presence 
of a clear periodicity immediately suggests a binary nature 
for the system. However, there are stars which show more 
complicated RV curves, and the interpretation of their spec- 
tral variability is not straightforward. HD 188209 (09.5Iab) 
is one of those objects. Garmany et al. (1980) have con- 
cluded from three spectra that this star is probably not 
a binary, and that the RV variations must be attributed 
to atmospheric motions. This conclusion was supported by 
Musaev & Chentsov (1988). However, based on 21 measure- 



ments Stone (1982) has concluded that HD 188209 can be 
considered as a spectroscopic binary with a period 57 days 
and small semiamplitude. More recently, FuUerton, Gies & 
Bolton (1996) included HD 188209 in a large sample of stars 
investigated on the presence of line profile variability (LPV) 
and found LPVs only in He I 5876 A. However, they did not 
flag HD 188209 as a velocity variable (their Table 10). 

The binarity of many O supergiants has been proposed 
recently by Thaller (1997). The fact that binaries have a 
higher incidence and an Ha emission strength in post-MS 
stages may indicate that wind interactions are a common 
source of emission in massive stars. In other words, even in 
cases where RV measurements are not available, the presence 
of Ha emission in a spectrum could be linked with colliding 
winds. One needs to study orbital phase variations in the 
Hq profile in order to be sure that the latter is due to collid- 
ing winds instead of some other mechanism. Note that HD 
188209 is an X-ray source detected by ROSAT (Berghoefer, 
Schmitt & CassineUi 1996). 

In this paper we focus on the high-resolution spectro- 
scopic data of HD 188209. Our observations can possibly 
account for the small semi-amplitudes and eccentric orbits 
of this binary candidate since they have been accumulated 
at different periods over a long baseline. 
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2 OBSERVATIONS 

The observations have been carried out in different runs (Ta- 
ble 1) using the the Coude Echelle Spectrometer (Musaev 
1993) at the 1-m telescope of the Special Astrophysical Ob- 
servatory of the Russian Academy of Science. Most of the 
spectra have a signal-to-noise ratio S/N>100 per resolu- 
tion clement, and an average resolution R = 40000 in the 
wavelength region 4400-7000 A. Preliminary reduction of 
the echelle spectra CCD images was made using the dech 
code (Galazutdinov 1992), which allows the flat-field di- 
vision, bias/background subtraction, one-dimensional spec- 
trum extraction from two-dimensional images, excision of 
cosmic-ray features, spectrum addition, correction for dif- 
fuse light, etc. Numerous bias, flat-fleld have been obtained 
every night. Each image was subject to a bias-frame subtrac- 
tion and fiat-field division using nightly means. Comparison 
exposures of a Th-Ar lamp were taken for each stellar spec- 
trum. The control measurements of interstellar Can and 
Nai D lines revealed a small scatter of 0.8 km s~^ (1 a). 
However, the 1 a dispersion of the velocity of the DIB was 
1.5 km s~^. These interstellar lines have been used to align 
all the spectra in the time series accurately. As an indicator 
of the overall precision of our measurements we have adopted 
1 a dispersion 1.5 km s~^ of the velocity of the DIB. All 
stellar absorption linos exhibited variations about their re- 
spective mean velocities at least 2-3 times the dispersion of 
the DIB velocities. The mean rms obtained from different 
dispersion curves was at least 0.003 A. The Coude Echelle 
Spectrograph was not a subject to mechanical and/or ther- 
mal instabilities. 

We used a 580x530 (pixel size 24 x 18 lira) CCD cam- 
era in all runs except 8 and 9. The last run was carried out 
using the Coude Echelle Spectrograph (Musaev, 1999) at the 
2-m telescope located in Terskol (North Caucasus, Russia). 
The CCD used in last two runs had a larger matrix (WI 
1242x1152 pixel with pixel size 22.5 x 22.5 /xm) allowing a 
coverage in a single exposure of the region ~ 3500-10100 A 
with almost the same resolution. 



3 PHOTOMETRY 

The photometry of HD 188209 obtained with Hipparcos 
(ESA 1997) is presented in Figure 1. The approximate re- 
sponse curve for the Hpdc passband (see van Leeuwen et al. 
1997) is very extended with a maximum at ~ 4500 A. Our 
observations did not reveal a strong variability in the equiv- 
alent widths of the strongest lines in the spectrum of HD 
188209, which suggests that the variations detected by Hip- 
parcos in the Hpdc band are due to changes in the continuum 
flux. The mean value of Hpdc=5.605 has a standard devia- 
tion 0.0155. The level of photometric variability is significant 
and cannot be ascribed to standard errors in the Hpdc mag- 
nitude of the order of 0.005 (van Leeuwen et al. 1997). Tests 
on the periodicity of the photometric variations were per- 
formed but no convincing period has been found. However, 
considering the large gaps in different Hipparcos measure- 
ments we cannot definitely rule out short-term photometric 
periodicity. 



Table 1. Journal of Observations. The second column gives the 
universal plus exposure time and in third column is the spectral 
region. The S/N ratio is given at 6600 A. 



Run 


Date 


UT + EXP 


Region 


S/N 




iHH Tnm wl 






(A) 




1 


14.06.93 


20h 


35m-|-30m 


4800-7500 


100 




r^o 

07.07.93 


20h 


40m-|-30m 


442()-o71U 


100 




f\Ci r\o 

08.07.93 


19h 


25m-|-20m 


4350-4510 


200 










6540-6720 




2 


04.09.93 


21h 


53m-|-30m 


4350-4510 


200 










6540-6720 






C\f\ 

05.09.93 


23h 


55mH-20m 


4350-4510 


200 










6540-6720 






C\f^ C\f\ 

06.09.93 


20h 


40m-|-20m 


A^^ C\ A^ ^ r\ 

4350-4510 


200 










6540—6720 






c\f\ 

07.09.93 


20h 


07m-f45m 


4350-4510 


200 










6540-6720 






nr\ r\€\ r»o 

09.09.93 


22h 


40m-t-45m 


4350—4510 


200 










6540-6720 




3 


07.1U.93 


17h 


46m-(-60m 


4370-6720 


150 




08.10.93 


17h 


35m-|-60m 


4370-6720 


150 




09.10.93 


17h 


40mH-60m 


4370-6720 


150 




10.10.93 


18h 


32m-t-60m 


4370-6720 


150 


4 


31.10.93 


16h 


22m-|-30m 


4330-4520 


200 










6540-6720 






01.11.93 


20h 


32m-|-30m 


4330-4520 


200 










6540-6720 






r\ o -1-1 /"\ o 

02.11.93 


16h 


18m+45m 


A A c\r\ r'C'jr\c\ 

4400-6800 


100 


5 


24.11.93 


17h 


56m+45m 


4420-7000 


100 




25.11.93 


17h 


02m-|-60m 


4420—7000 


150 




27.11.93 


19h 


15m-|-60m 


4420-7000 


150 




28.11.93 


16h 


50m-|-60m 


4420—7000 


150 




29.11.93 


16h 


41m-|-60m 


4420-7000 


150 




30.11.93 


15h 


50m-(-60m 


4420—7000 


150 




02.12.93 


16h 


22m+60m 


4420-7000 


150 




05.12.93 


17h 


08m+60m 


4420-7000 


150 


6 


16.09.94 


17h 


27m+60m 


4420-7000 


150 




"1 r\c\ c\ A 

17.09.94 


17h 


26m-|-60m 


4420—7000 


150 




18.09.94 


17h 


18m-|-60m 


4420—7000 


150 




19.09.94 


17h 


36m-f 60m 


4420—7000 


150 




20.09.94 


21h 


21m-|-60m 


4420-7000 


150 




21.09.94 


18h 


01m-|-60m 


4420—7000 


150 


7 


13.10.94 


17h 


00m-|-60m 


4350-6710 


150 




15.10.94 


16h 


40m+60m 


4350-6710 


150 




16.10.94 


20h 


20m+60m 


4350-6710 


150 




17.10.94 


18h 


01m+60m 


4350-6710 


150 




25.10.94 


18h 


00m-(-60m 


4350-6710 


150 




27.10.94 


19h 


06m-|-40m 


4350-6710 


100 




29.10.94 


16h 


36m-|-40m 


4350-6710 


100 


8 


15.05.97 


OOh 


35m-f45m 


3380-10060 


200 




15.05.97 


21h 


56m-|-60m 


3380-10060 


250 




16.05.97 


21h 


09m+60m 


3380-10060 


250 




30.08.97 


19h 


03m+30m 


4420-6800 


100 


9 


08.03.98 


02h40m-|-30m 


3560 10060 


200 




09.03.98 


01h50m+45m 


3560-10060 


250 




10.03.98 


02h02m+45m 


3560-10060 


250 




11.03.98 


01h32m-|-20m 


3560-10060 


200 
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Figure 1. Hipparcos photometry of HD 188 209. 



4 FUNDAMENTAL PARAMETERS OF 
HD 188209 

4.1 Plane— parallel analysis 

To analyse the spectrum of HD 188209 we first determine 
the rotational velocity from the width of 11 metal lines of 
C, N, O, Si, Mg and Ca, adopting a Gaussian instrumental 
profile with a FWHM of 0.13 A. The resulting value was a 
projected rotational velocity of 82.0 ± 8.5 km s~^, in good 
agreement with the value of 87 km s~^ reported by Penny 
(1996), and between those given by Conti and Ebbets (1977) 
and Howarth et al. (1997) who give 70 and 92 km s~^, re- 
spectively. 

The method followed in determining the stellar param- 
eters from the spectrum using NLTE, plane-parallel hydro- 
static model atmospheres has been described in detail by 
Herrero et al. (1992, and references therein). Briefly, we de- 
termine, at a fixed helium abundance, the gravity that best 
fits the different H and He profiles at a given temperature 
for a set of temperatures. If the abundance is right, the lines 
in the Toff-log g diagram will ideally cross at a point, giving 
the stellar Toh and log 3. Usually, they form an intersection 
region, whose central point is taken as giving the stellar pa- 
rameters, and whose limits give the adopted error. If the 
lines do not cross at any point, the helium abundance is 
changed. The helium abundance giving the smaller intersec- 
tion region for all profiles is the one selected. The center 



of the intersection region is taken again as that giving the 

stellar parameters. 

Recently, McErlea n, Le nnon & Dufton ( 199^ ) and 
Smith & Howarth (1998) have shown that different 
He I lines give different helium abundances in the region 
of the Teff- log g diagram occupied by HD 188209. They at- 
tribute this to the effect of microturbulence and show that a 
value of around 10 km s~^ is appropriate for bringing most of 
the He I lines into agreement. Thus, we adopt this value for 
HD 188209 and carry out the analysis in the way described 
above. 

With this method we have determined the stellar pa- 
rameters of HD 188209. We obtained Toff^ 31 500 K ±loo° , 
log (7= 3.0 ±0.1 (uncorrected for centrifugal force; because 
of the low rotational velocity we neglect this small correc- 
tion here) and e= 0.12 ±0.03 (the abundance of helium with 
respect to the total abundance of hydrogen plus helium, by 
number; the solar abundance is e— 0.09). The final fits for 
the H and He lines are shown in Figure 2. 

The larger errors towards higher temperatures is due 
to a small difference in the two spectrograms available for 
Hell 4541 A. Using the second one, we would have obtained 
a Teff of 32 000 K, all other parameters remaining the same. 
For this reason, we have enlarged the error bar in this direc- 
tion. Remember that the errors given are formal errors, in 
the sense that they express the uncertainty in the fit using 
the models described above. 



As has been show by Herrero (1994) the metal opacity 
in the UV (line blocking) could also affect the value of the 
parameters determined. However, at the relatively low tem- 
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Figure 2. The fit to the HD 188209 lines using unified models, a) H7; b) He I 4471; c) Hei 4387; d) He I 4922; e) Hell 4200, and f) 
Hen 4541 A (see text for details). 



The spectral variations of the 0-type runaway supergiant HD 188209 

b) 





0.60 



4330 4340 4350 

WAVELENGTH IN ANGSTROMS 



4360 



4465 4470 4475 

WAVELENGTH IN ANGSTROMS 



4480 






0.80 

4190 



0.70 



1.05 



1 .00 



0.95 - 



4915 4920 4925 4930 

WAVELENGTH IN ANGSTROMS 



> 
< 



0.90 - 



0.85 



0.80 




4195 4200 4205 

WAVELENGTH IN ANGSTROMS 



4210 



4535 4540 4545 4550 

WAVELENGTH IN ANGSTROMS 



6 G. Israelian et al. 



perature of HD 188209 the effect would be minor, moving 
Teff towards higher temperatures within the error box. 

With the stellar parameters given above we can deter- 
mine the radius, luminosity and mass of HD 188209 as de- 
scribed in Herrero et al. (1992). For A/v= —6.0 mag (Mboi 
= —9.0 mag) given by Howarth & Prinja (1989) we ob- 
tain R/Rq= 20.9, log(L/L0)= 5.59 and M/Mq= 16.6. 
The errors are again as in Herrero et al. (1992): ±0.06 in 
log(i?/i?0), ±0.16 in log(L/LQ) and ±0.22 in log(M/MQ). 

HD 188209 does not formally show the helium discrep- 
ancy, as the solar helium abundance is within the error bars. 
However, it shows the mass discrepancy: the mass derived 
from the plane-parallel spectroscopic analysis is, even in- 
cluding the error bars, much lower than the o ne der ived 
from the evolutionary tracks from Schaller et al. (1992). 



4.2 Unified model analysis 

After having the parameters from the plane-parallel anal- 
ysis, we can try to use a spherical, non-hydrostatic model 
atmosphere in order to improve the already derived param- 
eters and also to obtain the mass-loss rate. In a supergiant 
like HD 188209 this can have an important impact on the 
final parameters. Usually, it is also assumed that this will 
contribute to the reduction in the mass discrepancy. 

The unified code we use is that recently developed by 
Santolaya-Rey, Puis & Herrero (1997). The reader will find 
all the details therein, but for our present purposes we men- 
tion that the code uses spherical geometry, with a /3-velocity 
field, and treats the wind and the photosphere in a unified 
way. It also makes use of the NLTE-Hopf functions. Stark 
broadening is included in the formal solution, and the model 
atoms are the same as in the plane-parallel case (slightly 
adapted for the new program). 

We begin by estimating the mass-loss rate from the 
Ha profile, and then, with this meiss-loss rate, we try again 
to find the best gravity (from the H7 wings), effective tem- 
perature and He abundance (from the He ionization equi- 
librium). With the new parameters, we again try to fit 
the Ha by varying the mass-loss rate, and so on. In the 
whole process we take the wind terminal velocity from Haser 
(1995), who gives 1700 km s"^ Note that Howarth et al. 
(1997), give a similar value of 1650 km s~^. 

We have used for the analysis the same lines as in Fig. 2, 
but have included He I 4471 A instead of H/3(that improved 
as does H7) to show the dilution effect (see below) . In Fig. 3 
we show the fit of the H7, He I and Hell lines. The adopted 
parameters are now TaB= 31 500 K, log (7= 3.00 and e— 0.12, 
i.e. the same parameters as in the plane-parallel case, even 
for the gravity (again adopting a microturbulent velocity of 
10 km s~^). Thus, the unified models do not contribute in 
this case to changes in the mass discrepancy found with the 
plane-parallel models (nor in the He abundance). We can 
also see in Fig. 3 that the He l 4471 line show s the well known 
dilution effect mentioned by Voels et al. (198!:), although 
the other He lines fit perfectly well. This effect merits an 
explanation. 

The fitting of the Ha profile needed to derive the mass- 
loss rate cannot be done properly. The profile is highly vari- 
able, and we have adopted a qualitative approach: we have 
simply tried to give upper and lower limits to the mass-loss 




0.851 L 

6550 



Figure 4. The Ha profiles of HD 188209 observed on 1993 
September 7 enclosed by difli'erent theoretical profiles calculated 
for log M= -5.70, -5.80 and -5.90 (dotted, dashed and dash- 
dotted lines, respectively) for the parameters given in text. The 
selected observed Ha profile is neither an average nor represen- 
tative, and the figure merely illustrates the approach we have 
followed. 



rate. As an example, we illustrate the procedure in Fig. 4, 
where we show one of the profiles with various mass-loss 
rates for the stellar parameters given above. The theoreti- 
cal mass-loss rate values shown in Fig. 4 are log M= —5.70, 
—5.80 and —5.90 (with M in Mp^/yr ). Thi s situation is sim- 
ilar to that found by Herrero et al. ( ]199^ ) for Cygnus X-1. 
The profile shown cannot be adopted as either an average 
or a representative one, as the profile varies a lot. The fig- 
ure is only for illustrative purposes. An average logarithmic 
mass-loss rate for HD 188209 would be between —6.0 and 
—5.7. These values are also in agreement with the profiles in 
Fig. 3, where the adopted mass-loss rate was log M=— 5.80 
(corresponding to 1.58 10~^ Mq/jt). We should point out 
that other values of the mass-loss rate give worse fits to the 
H7 and He I lines, by strongly modifying the line cores. 



5 RADIAL VELOCITY VARIATIONS OF THE 

ABSORPTION LINES 

The radial velocity variations in stars can have instrumen- 
tal, internal (atmospheric) and/or external (Keplerian) ori- 
gin. Instrumental effects in our measurements are minimized 
due to the high resolution and high S/N of the data pre- 
sented. The internal accuracy achieved for the wavelength 
calibrations is of the order 1.5 km s~^ as derived from the 
scatter of measured radial velocities of interstellar and tel- 
luric lines in the spectra. Note that all former studies of 
HD 188209 (except for five spectra obtained by FuUerton 
et al. 1996) were based on the photographic spectra. Atmo- 
spheric pulsations of early-type stars have been the subject 
of extensive studies (e.g. Burki 1978; Bohannan & Garmany 
1978, Kaufer et al. 1997, FuUerton et al. 1996) and have 
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Figure 5. Mean radial velocities of different groups of lines ver- 
sus TEE. A linear fit (0.1527,-17.05) provides a non-paxametric 
correlation coefficient 0.73. 
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Figure 6. Standard deviations of mean radial velocities of differ- 
ent groups of lines versus TEE. A linear fit (—0.049,8.474) pro- 
vides a non-parametric correlation coefficient 0.76. 



been reviewed by Baade (1988, 1998). Gies (1987) has com- 
pared the velocity distribution and binary frequency among 
195 Galactic O-typc stars (cluster and association, field and 
runaway) and found a deficiency of spectroscopic binaries 
among field stars (and especially among the runaway). The 
most comprehensive radial velocity studies of OB and BA 
supergiants to date are those of FYiUerton et al. (1996) and 
Kaufer et al. (1996; 1997). Radial pulsation periods with 
P<5'* have been predicted for O-type supergiants (Burki 
1978; de Jager 1980; Levy et al. 1984) . The presence of pul- 



sations or random motions in stellar atmospheres results 
in complex velocity curves for different spectral lines due 
to stratification effects (Abt 1957). This is in contrast to 
\ Keplerian motions where all the lines vary synchronously 

- with time (Ebbets 1979; Garmany et al. 1980). However, 
in many cases the amplitude of the RV variations docs not 

- exceed 25-30 km s"'^, and it is very difficult to distinguish 
J whether these variations are of an internal or an external na- 

- ture. Additional difficulty comes from the possible presence 
of non-radial pulsations (NRP). Variable profiles (LPVs) 

- have been detected in many narrow-line supergiants (e.g. 

- Baade 1988; Kaufer et al. 1997; FuUerton et al. 1996) and in 

- some cases they correspond to the radial velocity variations 

- measured in photographic spectra. However, in contrast to 
broad-line supergiants, clear evidence of NRPs in narrow- 

- and intermediate-line O-type supergiants has not yet been 
found. Although LPVs are a common occurrence among the 
O-type stars, some of them (like HD 34656; Fullerton, Gies 

, Bolton 1991) show a variability which consists of cyclical 
fluctuations in radial velocity due to pulsations in a funda- 
mental mode. 

Careful inspection of our time series showed that all 
absorption lines varied in position over the course of a day. 
Wc have detected a few asymmetric profiles of He linos but 
there were no signatures of moving features. It is however 
possible that the S/N ratio of our data is not high enough to 
trace LPVs. Small distortions in profile shape axe typically 
less than 1% of the continuum strength and a minimum S/N 
ratio required would be at least 300 (Fullerton et al. 1996). 



5.1 The velocity-excitation relationship 

We have selected unblended lines by utilizing theoretical 
synthetic spectra computed for the model atmosphere of 

- HD 188209 and measured RVs by fitting a Gaussian to 

- the line profile. The following groups of lines have been se- 

- lected: Hei (4921, 4471, 4713, 5015, 5047, 5876, 6678 A), 

- Hell (5411, 4541, 4686 A), Siiv (4654, 4631 A), Cm (4650, 
_ 4647, 5696 A), Sim (4567, 4552, 4574 A), Civ (5801, 
_ 5811 A), Niii (4514, 4510, 4518, 4523 A), On 4661 A and 

Mgll 4481 A. The average RVs computed for each of these 
groups of lines are listed in Table 2. It is normally £is- 

sumed that lines of different total excitation energy (TEE, 

! Qoionization energy plus excitation energy of the lower level) 
form in different layers of the atmosphere (Hutcliings 1976). 
However, it is not clear whether the stratification exists in 
a dynamically active, pulsating atmosphere. One can as- 
sume that the time scale of dynamical processes (pulsations, 
stochastic motions etc.) is much less than the time necessary 
for the establishment of radiative equilibrium. Thus, while 
pulsating, the atmosphere is supposed to pass through a 
chain of hydrostatic stages. The assumption that TEE cor- 
relates with line formation depths can be tested as well. 
Wc applied our plane parallel models to compute formation 
depths of the line cores of the He I and He ll lines. It appeared 
that the core of the strong Hell 4686 A line forms much 



closer to the surface (at column mass ~ 0.01 gcm~ ) than 
any of the He I lines. However, this was not the case with 
the two other Hen lines (5411 and 4541 A) wc used in our 
study. Similar tests have been carried out for other groups 
of lines but assuming an LTE line formation. These exer- 
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Figure 7. Lomb-Scargle (top panel) and CLEAN periodograms of tlie average RV for each date obtained by averaging tiie RVs of all 
the groups. 



ciscs suggest that before combining lines in difTcrcnt groups 
and computing their average RVs, one has to be sure (of 
course under the assumption that our plane-parallel models 
are applicable) that their depths of formation are similar. 

The velocity-excitation relationship found for many hot 
supergiants (Hutchings 1976) exists also in HD 188209. In 
Figs. 5 and 6 we present plots of mean RV versus TEE and 
standard deviations of the mean RV versus TEE for the dif- 
ferent groups of lines computed for all dates (the last two 
lines in Table 2). These plots exclude a pure Keplerian mo- 
tion as the only cause of the RV variations. Pulsations and 
stochastic motions (intrinsic wind variations caused by some 
hydrodynamical instabilities) in the wind can bring to the 
RV variations as well. If we suppose that stochastic varia- 
tions are not important, then the existence of a standard 
deviation-TEE relationship would suggest that deeper lay- 



ers in the atmosphere pulsate with smaller amplitudes. The 
amplitude of pulsations increases when approaching to the 
surface. 



5.2 The periodicity of the radial velocity 

variations 

The period search was carried out with help of the PERIOD 
package (Dhillon & Privett 1997) of the starlink software. 
The following strategy was applied when looking for a pe- 
riodic signal in the RV curves of different groups of lines. 
Due to the large gaps (especially between runs 5 and 6) in 
our observations, we first decided to study each of the runs 
5, 6 and 7 separately. The clean algorithm (Roberts et al. 
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Table 2. Radial velocities of different groups of lines. 



Date He I Heii Siiv Cm Sim Civ Nm On 4661 A Mgii 4481 A 

TEE (eV) 21 75 57 54 35 85 57 36 8.5 
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1987) was employed to cover a space of loop gains from 0.2 
to 0.6 and the number of iterations from 10 to few hundreds. 
The convergence of the periodograms was achieved for the 
majority of groups of lines in all three runs. The mean fre- 
quency suggested by most of the groups in all three runs is 
0.44±0.05 days"'^ (2.24 days). However, this period is very 
close to the Nyquist frequency (1/(2 x Smallest Data Inter- 
val)) of the data and might be misleading. 

We have also looked for periodic signals in the combined 
data of all groups obtained in all runs (Table 2). The max- 
imum and minimum frequencies were set to 100 and 0, re- 
spectively. A CLEAN analysis of the time series of the major- 
ity of groups revealed a frequency of 0.51±0.1 days^^ (1.95 
days). The gain factor was 0.1 at the first iteration, then 
was decreased by 15-20 iteration until stabilization. The av- 
erage RV for each date obtained by averaging the RVs of 
all the groups revealed a frequency of 0.47±0.12 days~^ (2.1 
days). Again, both frequencies arc very close to the Nyquist 
frequency and we should discard them. We must point out 
that our periodograms did not show any peaks at frequen- 



cies smaller than 0.4 days~^. The next strongest peak which 
appeared in our periodograms was near 0.156±0.15 days~^ 
(6.4 days). Clearly this period is not affected by sampling. 
We have also analysed the RV data using the Lomb-Scargle 
method (Lomb 1976, Scargle 1982) which allows to compute 
statistical probability of peaks in periodograms. To ensure 
reliable significance values, the minimum number of permu- 
tations was set 100. The probability that the period is not 
equal to 6.4 days was always less than 30 %. The peak at 6.4 
days appears in all periodograms but given its significance 
value, we cannot definitely rule out its non-physical nature. 
In Figs. 7 and 8 we show the clean and the lomb-scargle 
periodograms and the fitting of a sin curve to folded data, 
respectively. 
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Figure 9. The variety of Ha profiles observed in difi'erent runs. A vertical dashed line indicates zero velocity. 



6 VARIABILITY OF Ha 

All hot supergiants have variable Ha profiles in their spec- 
tra (Rosendahl 1973) . The shape of the Ha may vary from 
P Cyg to inverse P Cyg, double-peaked, pure absorption 
and/or emission (Ebbets 1982) with typical time-scales of 
the order of days. The nature of this variability is not 
yet understood. The existence of variable asymmetric out- 
flows/infalls of matter and some corotating structures re- 
lated to surface inhomogonitics and possible magnetic fields 
have been proposed for BA-type (Kaufer et al. 1996) and O- 
type (Fullerton et al. 1996; Kaper et al. 1997) supergiants. 
In addition, there have been detailed studies of the rotat- 
ing giant loop in /3 Orionis (Israelian, Chentsov & Musaev 
1997) and the corotating spiral structures in HD 64760 and 
HD 93521 (Howarth et al. 1998; Fullerton et al. 1997). It is 
of course very difficult to distinguish binary systems from 



single stars without understanding the nature of Ha vari- 
ability. As Thaller (1997) suggests, the Ha can suffer some 
pecuhar variability due to the colliding winds in a binary 
system. 

The time evolution of Ha profiles in three different runs 
is shown in Figure 9. The average Ha profile consists of three 
components, a central emission accompanied by blue and red 
absorptions. We have not observed a single Ha profile with- 
out a central reversal. The emission is not always centered 
exactly on the rest wavelength but is varying. It may ap- 
proach the continuum level, go above it and decrease rapidly 
in strength. Apparently the time-scale of the Ha variability 
is at least one day. The 5*'' run has been divided into two 
parts (runs 5a & 5b) with four successive nights in each. 
The Ha variability is observed over a wide range from about 
-400 to 200 km s~^ 

We have already seen in Section 4.2 that our spherical 
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Figure 10. The radial velocities of blue and red absorption components versus residual intensity of the central emission. 



unified models can account for the central reversal (or at 
least set upper and lower limits of the mass- loss). Thus, we 
know that the central emission forms in the expanding en- 
velope and accounts for the filling-in efi'ect observed in Ha, 
H/3 and some other lines. The filling-in effect observed in 
H/3 correlated perfectly with the strength of the Ha central 
reversal. The Ha wings originate deep in the atmosphere, 
whereas the central reversal comes from the thin layers of 
the envelope. One cannot use the term "underlying pho- 
tospheric absorption line" since the central reversal is not 
emitted by a detached layer far from the photosphere. It is 
important to stress that we deal with a single line formed 
in a unified model atmosphere. A variable amount of in- 
cipient emission can be due to density (or radius) variations 
in the outer atmosphere. However, these variations are not 
expected to produce an asymmetry as long as we deal with 
spherically symmetric mass-loss. Our observations indicate 



that the velocity of the central emission varies as well. We 
do expect RV variations in the central emission (the upper 
atmosphere) since we know that the RVs of all absorption 
lines vary. The amplitude of the RV variations of Mg ll 4481 
A can reach 30 km s"^ (Table 2) and 40 km s"^ in H/3. 
Thus, it is not unusual for the amplitude of RV variations 
of the Ha central emission to reach 50 km s~^. A period 
analysis of the RV curves of the central emission resulted in 
the detection of quasi-periodic variability with a frequency 
0.42±0.14 d"^ (2.35 days). It turns out that the RV curves 
of all groups of lines vary in phase. However, the RV curve 
of the Ha was shifted half phase relative to all groups. 

We have measured the RVs of the blue and red absorp- 
tion components of the Ha and plotted them against the 
residual intensity of the central emission reversal. This plot 
(Fig. 10) shows a correlation with a large scatter due to the 
RV variations of the central reversal. This kind of correla- 
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Figure 11. The phase diagram of the EWs of Ha. 
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Figure 12. Dynamical phase spectrum of Ha. 

tion can be expected when the central reversal is moving up 
and down relative to a local continuum. We found similar 
correlations between the RV of the central emission and the 
residual intensities of the red and blue absorptions. Appar- 
ently all the changes observed in the Ha wings at velocities 
V < — 100 km s^^ and v > +100 km s^^ are due to the varia- 
tions of the central reversal. We do not anticipate such large 
RV variations deep in the atmosphere where these wings are 
formed. The conclusion is that the overall shape of the Ha is 
determined by the central emission. 

Wc have performed a period search of the integrated 
equivalent width (EW) of the Ha data set and found a 
maximum in power at frequency 0.22 day"'^. Figures 11 and 
12 show the phase diagram for the period 4.41 days and 
a grey-scale representation of the phase spectrum, respec- 
tively. A phase spectrum represents a two-dimensional case 
of the CLEAN algorithm where each velocity bin of the Ha is 
treated as a time series of the Ha intensity. 



7 DISCUSSION 

Our target belongs to the group of stars for which the ex- 
istence of a compact companion has been proposed in the 
literature. The tcisk of disproving or confirming the binary 
nature of the system can be tackled only if sufficiently accu- 



rate analysed observational data are available. In this paper 
we used state-of-the-art models of atmospheres to determine 
the fundamental parameters of HD 188209. 

To establish the presence of a possible companion we 
have studied the RVs of absorption lines by combining them 
in different groups. Fourier analysis based on the iterative 
CLEAN algorithm wEis used to search for periodic variability. 
Unfortunately the time coverage of runs 1-4 and 8-9 was too 
sparse to set constraints on their time-dependent behaviour. 
For this reason wc first analysed a few runs separately and 
then utilized the clean algorithm to search for periods in 
a whole data set. The highest peak in the Fourier power 
spectrum was centered near the frequency 0.156 day~^ (6.4 
days). 

The 6.4 days period can bo due to the binary nature 
of the system if one assumes very small and unlikely values 
for the mass ratio (q<0.1). Taking the values derived in this 
article {M/Mq= 16.6, R/Rq=20.9) and assuming q=0.1, we 
obtain for the Roche radius and for the major semi-axis of 
the binary orbit 15 -R0 and 25 -R0, respectively. This simple 
estimate shows that an O supergiant can hardly fit within 
the orbit because its Roche radius would be less than the 
stellar radius. Even if it would fit, the tides in such a tight 
binary would be very strong making the star to speed up 
quickly until the rotation period matches the orbit. Even if 
we assume that the system is very young, it's hard to explain 
that the putative orbital period is 2 times shorter than the 
rotational period (13 days). The second difficulty with the 
binary interpretation comes from the variability of Ha. A 
TVS (temporal variance spectrum, showing the extent and 
distribution of statistically significant profile variability) has 
been computed recently (Baade 1998b, Kaper et al. 1998) 
for 15 spectroscopic binaries and it was found that all they 
show a characteristic double-peaked profile. This is due to 
two Ha absorption/emission profiles moving in a composite 
spectra. In our case the Ha profile is splitted because of 
the central emission coming from the lower wind. The last 
argument comes from the clear relations between excitation 
energies and radial velocity amplitude and excitation energy 
and mean radial velocities and from the model atmosphere 
atmosphere calculations. The latter is a good discriminant 
between internal variations (pulsations & wind instabilities) 
and Keplerian motions. In a binary system one would expect 
all lines to have the same amplitude independent on their 
TEE. 

It has been known for a long time (Abt 1957) that the 
quasi-pcriodicity in hot supcrgiants might be ascribed to 
radial pulsations. A simple relation (Burki 1978; de Jager 
1980) can be used to estimate the period of radial pulsation, 

logPfund = 10.93-0.5 log(M/M0)-O.38Mboi-31ogTefi(l) 

Using the values of parameters obtained in Section 4 we ar- 
rive at Pfund=l-75 d. Note that the form of the relation (1) 
depends on the stellar evolutionary models and the input pa- 
rameters; both are subject to large errors. In particular, note 
that we found no large differences in the parameters deter- 
mined with plane-parallel and unified model atmospheres. 
Nevertheless, Levy et al. (1984) have pointed out that pe- 
riods a factor of 1.5 longer than the corresponding periods 
of the radial pulsations can be ascribed to non-radial pul- 
sations. This means that a factor of two difference between 
the evolutionary and the spectroscopic masses can easily re- 
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suit in the mis-identification of the pulsating mode. Another 
difficulty has been pointed by the referee of the article. A 
more sophisticated approach shows (Unno et al. 1979) that 
f-mode pulsation (which is the lowest-frequency mode sup- 
ported by radial pulsation) periods are about 10 times larger 
than the one suggested by a period-luminosity relation. In 
any case, the theoretical period of 1.75 days is very close 
to the Nyquist frequency of our data which means that we 
have a little chance to identify it in our data set even if it 
exists. 

Our data not allow to distinguish between pulsations 
and stochastic variations of the stellar wind. It is also quite 
possible that we have a combination of both effects. 

Note that the projected rotational period of this star 
(~ 13 d) is much longer than any of the quasi-periods found 
in this paper (but of course our runs do not cover a whole 
rotation cycle). The surface features (if any) will always be 
visible on the projected disc of the star independently of 
the inclination angle. Thus, any periods due to the rotation 
of these features must correspond directly to the rotation 
period. We do not find any peaks in the power spectra at 
~13 d and this leads us to discard rotational modulation as 
a possible explanation of the RV variations reported here. 

The quality and the sampling of our data do not allow a 
careful study of the line asymmetries, moving components (if 
later exists) and/or long-term spectroscopic variability to be 
made. It is quite possible that the non-sinusoidal character 
of the RV curve for 6.4 days period (Fig 8) is caused by 
some disturbances due to the NRPs and/or moving features 
in the profiles plus any stochcistic instabilities of the wind. 
New monitoring with much higher S/N may allow NRPs, 
multimode pulsations and clearly separate a sinusoidal curve 
of the radial pulsations to be revealed. However, we found 
convincing evidence that the atmospheric motions cannot 
be ascribed solely to Keplerian motions and probably are 
not of a binary origin. 
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